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Asymptotic behavior of the Tomonaga-Luttinger liquid in the vicinity of the phase-separated region 
is investigated in the one-dimensional t-J model, to study the universal property of the c = 1 
conformal field theory with U(l) symmetry near the K —* oo instability. On the analogy of the 
spinless fermion, we discuss that the compressibility behaves as k oc [Jc — J)~^ , and that the Drude 
weight is constant and changes to zero discontinuously at the phase boundary. This speculation is 
confirmed by analyzing the finite size effect from the result of the exact diagonalization. 
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I. INTRODUCTION 

Phase separation in one-dimensional (ID) electron sys- 
tems has attracted attention of condensed matter physi- 
cists, since it is observed that the superconducting cor- 
relation is enhanced near the phase-separated |jjegion. 
Many electron systems show phase separationld'Bfla, of 
which the t-J model is the simplest. The Hamiltonian of 
the ID t-J model is written as 

n = -i^(4^c,+i^ + 4-HiaC*^) 

+J^{S^-Si+i-nm+i/4:), (1) 

i 

in the subspace without double occupancy. In general, 
unless some instability occurs, ID electrorusj/stems be- 
long to Tomonaga-Luttinger liquids (TLL)la'El, which is 
characterized by gapless charge and spin excitations and 
power-law decay of correlation functions. The charge 
part of the excitation is described by the c = 1 confor- 
mal field theory (CFT) with U(l) symmetry, while the 
spin part is described by the c = 1 CFT with SU(2) sym- 
metry. The dominant correlations are determined by a 
single exponent Kp. In the phase diagram of the ID t-J 
model obtained by Ogata et al.u, the charge or the spin 
density fluctuations are dominant {Kp < 1) for the small 
J/t region, while J/t is increased, the superconducting 
fluctuations become dominant {Kp > 1). Finally, the 
compressibility diverges at Jc/t — 2.5-3.5 {Kp -^ oo) 
and the system goes into the phase-separated state for 

J > Jc. 

As a typical instability of the c = 1 U(l) rCFT, 
Berezinskii-Kosterlitz-Thouless (BKT) transitionQ has 
been well known. In electtpn systems, the Mott tran- 
sition belongs to this typeB. However, the asymptotic 
behavior of the c = 1 U(l) CFT with K ^ oo has not 
been fully investigated. In ID spin systems, this type 
of instability corresponds to the transition between the 
massless XY phase and the ferromagnetic phase. In this 
article, we argue the asymptotic behavior in the vicinity 



of the instability K -^ oo, relating it with the ID spinless 
fermion system. Then we apply our speculation to the 
ID t-J model by analyzing the finite size effect from the 
result of the exact diagonalization. 

This paper is organized as follows. In section y|, we 
discuss an instability of the c = 1 CFT with U(l) sym- 
metry, and review several physical examples. In section 
III, we numerically calculate physical quantities in the ID 
t-J model to examine the argument. Finally, a summary 
is given in section 



IV 



II. PHASE SEPARATION AND THE C = 1 U(l) 
CFT 

The c — 1 CFT with U(l) symmetry is described by 
the Gaussian model 



27r 



dx 



2 V f dcj) 



^^(-n) +ifU^ 



(2) 



where 11 is the momentum density conjugate to 4>, 
[(l){x),'n{x')] ~ i5{x — x'), K is the Gaussian coupling, 
and V is the sound velocity. Its dual field 6 is de- 
fined as dx9{x) = 7rn(a;). We make the identification 
(j)=(j) + 27r/V2, e = e + 27r/%/2, which means the U(l) 
symmetry for the 9 field. The scaling dimensions for the 
operators 



are 



On,m = exp(mv 20) exp(ir7iv 20) 



_ 1 / 27^ ™ 



(3) 



(4) 



Therefore, the scaling dimension xo,m decreases to for 
K —t oo, which implies phase separation in electron sys- 
tems, or a ferromagnetic long-range order in spin sys- 
tems. 



For a finite L size system with periodic boundary con- 
ditions, thftnexcitation energies are related to the scaling 
dimensions^ 



AF (T\-—T 



(5) 



and the correction to the gpaund state energy is described 
by the conformal charge cE£l 



Eo{L) = eoL- — c, 
where eo is the bulk energy density. 



(6) 



A. Spinless fermion systems 

As the simplest physical model, we consider the spin- 
less fermion with nearest neighbor interactions {t-V 
model). 



^ = -^E(cI^«+i + 4 



iL-jy 



V 2_] niTiiJ, 



(7) 



For simplicity, we replace the parameters as A = F/2|<|. 
Then the metallic phase, which is described by the c = 1 
U(p-|CFT, becomes unstable and phase separation {K -^ 
cx))Ell occurs for the region of A > 1 at any fillings. 

The quantjiuji numbers in (0) correspond as (n, m) = 
(AD, A7V)y'yta, where AD denotes the number of par- 
ticles moved from the left Fermi point to the right one, 
and AiV is the change of the total number of particles. 
The selection rule under periodic boundary conditions is 



AiV 
AD = — - (modi). 



(8) 



Since the maximum values of the A7V and AD are of 
the order of the system size, and the excitation energies 
Ai?A_D,o = ttAD^uA'/L are limited by the band width 
(oc L), so that vK is expected to be constant near the 
critical point {K — )■ oo). 

It is thought that the transition between the metal- 
lic state and the phase-separated state is the first or- 
der, because the ground state of the phase-separated 
state (the state where the band is filled up or empty, 
i.e., AA^ = ±L/2 at half-filling), and the metallic state 
(AiV = 0) is thermodynamically different. Therefore, 
the excited energy corresponding to the phase-separated 
state from the metallic ground state should be 



ADo, 



±L/2 



TTV {L/2f 

T K 



oc (1- A)L, 



(9) 



which means v/ K oc 1 — A. 

In fact, using the Jordan- Wigner transformation, the 
Hamiltonian (0) at half-filling is mapped pOpto the S — 
1/2 XXZ model with the anisotropy AO. From the 
Bethe ansatz resultll3, we obtain 



l/K = - arccos A w -J2{\- /^), 

IT IT 



(10) 



and 



^^^.sin(arccos^^ 1^2(1^^. (11) 

2 TT - arccos A 2 ^ ^ ^ ^ ' 

Therefore, near A = 1, the sound velocity behaves v oc 
1/A', as expected. Note that the phase-separated state 
has the excited energy near A = 1 



ADo,±L/2 = |(l-A), 



(12) 



which is consistent with (||). Thus, near A = 1, the 
asymptotic behaviors are vK ~ Cansi.^v/K oc 1 — A, as 
expected. 



B. Electron systems 

In the case of ID electron systems, generally, the 
charge and the spin degrees of freedoms are separated, 
and the charge part is described by the c = 1 U(U-.GET, 



while the spin part is described by the SU(2) CF1 

^ 27rWc 2ttVs 

AD = Xr H x^. 

L L 

Here the scaling dimensions are given by 



(AiV,/2)2 



A' 



Kp{2ADc + AD, 



a;, = - I AiV, 



ANr, 



Ud^ 



(13) 

(14a) 
(14b) 



where the subscripts c, s denote the charge and the spin 
degrees of freedoms respectively, and AiV^ means the 
change of the number of the down spins. These quan- 



tum numbers are restricted 
periodic boundary conditional 



the selection rule under 



AD, 



AD, 



A7V, + AN, 



ANr 



(mod 1) 
(mod 1) 



(15a) 
(15b) 



Since both the spinless fermion (the S = 1/2 XXZ 
spin chain) and the charge part of the electron system, 
for example the ID t-J model, belong to the same uni- 
versality class, we expect that the asymptotic behavior in 
the limit of AT — > oo is the same for both cases. Therefore, 
the charge velocity and the exponent for the charge part 
are expected to behave as VcKp ~ Const., v^/ Kp oc Jc — J 
respectively. 

Next, we relate Vc and Kp to the compressibility 
and |tlje-|Drude weight. From the universal relations of 
TLLElEil, the compressibility is given by 



TT Vc 

2K, 



(16) 



so that the compressibihty is expected to behave like k cx 
{Jc — J)^^- On the other hand, the Drude weight is given 
by the relation 



D 



KpVc 



(17) 



therefore, it is constant (at least finite) at the phase 
boundary. In the phase-separated state, the system 
is separated into an electron rich and a hole rich re- 
gions, and the former is regarded as an antiferromagnetic 
Heisenberg chain with open boundary conditions, so that 
the dc conductivity in this non-uniform state is expected 
to be 0. This indicates that the Drude weight changes 
discontinuously at the phase boundary. 



C. Spin systems 

For a model of spin systems, we treat the XXZ spin 
chain 



/I — y \^ j ^-'7+1 ' ^ j ^ j ' ^Oj Oj^-^j 



(18) 



This model has a massless XY phase close to the A < 1, 
and a ferromagnetic phase for A > 1. 

As mentioned in the previous subsection, the Hamil- 
tonian ( p8| ) for S — 1/2 is obtained by the Jordan- 
Wigner transformation of the t-V model (|^). However, 
the boundary condition depends on the number of par- 
ticles in the spinless fermion system: S^_^_■^ = e^'^'^^ S^ , 
which changes the selection rule to the bosonic one, 

AA^ = integer, AD ~ integer. (19) 

Changing the selection rule, our argument for the spinless 
fermion also holds in this case. The fact that v oc 1/K 
explains the change of the dispersion curve from the type 
(jj K, v\q\ in the XY region to the lo (x q^ o\y the SU(2) 
ferromagnet, and that the ground state energy of the 
ferromagnet does not depend on the system size. 

For the general spin 5* case, in the XY phase the fully 
ferromagnetic states (5*1, — zLSL) have the excited en- 
ergy 

^Eo,±sL = S^L{1 - A)[l + 0(1 - A, l/L)], (20) 

(see Appendix A). Comparing this with the Gaussian 
model arguments (0), we obtain 



AE, 



0,±m 



X<' 



TTV m 



(21) 



Therefore, assuming v oc 1/K near A = 1, we obtain 
w, \/K oc vT^-^- And the energy gap between the fully 
ferromagnetic state and the one-spin flip state (5*1^ = 
±(5'L-l))is 



2SL- 1 



L 



-(1-A)«25(1-A), 



(22) 



which is consistent with the simple spin wave calculation. 
For the ferromagnetic region (A > 1), there is a large 
degeneracy 0{L) in the ground state under the special 
boundary condition, which reflects the iijwjariance of the 
Hamiltonian under the quantum groupEa. Physically, 
this degeneracy is related with the translational invari- 
ance of the domain wall. 



III. NUMERICAL CALCULATION 

To examine the above prediction, we diagonalizc the 
Hamiltonian of the t-J model (|l|) using the Lanczos 
method for 8,12,16,20 sites clusters at quarter flUing. In 
the finite size calculation, the compressibility k is given 
by 

L f EojL; N + 2) + Eo{L;N^2)^ 2Eo{L; N) \ ~^ 
""^IP \ 4 ) 

(23) 

where Eo{L;N) is the ground state energy of a system 
with size L and A^ electrons {n = N/L). We choose 
periodic boundary conditions for N = Am + 2 (m: in- 
teger) electrons, and antiperiodic boundary conditions 
for A^ = Am electrons. The reason of this choice is 
as follows. When the number of electrons are changed 
by 2 (AA^c — 2) keeping S^ — 0, the number of up 
spin changes by 1 (AA^g — ^)- Then, from the selection 
rule (15a), the possible value of ADc shifts by 1/2. The 
change of the boundary condition cancels this phase shift, 
and it makjea|-the ground state always singlet with zero 
moment umEju. 

The Drude weighttJo is given with the relation 



D 



L d^Eoi^) 



a$2 



(24) 



$=0 



where the flux-dependent ground state energy _Bo( 'I') is 
calculated by modifying the hopping part of (|l|): 

-t5]e^*/^4c,+i,+H.c.. (25) 



By fitting the energy difference i?o(*&) ~ £-0(0) as a func- 
tion A$2 _^_ ^^4^ ^^ g^^ ^^g^-^ A:$> \B\> 0. Therefore, 
the Drude weight and $^ is well approximated by the 
linear relation, and we can calculate the Drude weight 
by the difference of energies (e.g., at $ = and O.Itt). 
This is a natural consequence in the TLL, because ADc 
is modified as ADc + ^/2tt for $ 7^ 0, so that ( p^ ) 
gives only an 0($^) dependence if the ground state is 
singlet. The Drude weight is known exactly at J/t = 0|-as 
D = sin(7rn)/7r, and at J/t = 2 from the Bethe ansatzEa. 



We have checked the value at these points and confirmed 
the vaUdity of our calculation. 

As we expected, the inverse compressibility k~^ van- 
ishes linearly toward the phase boundary in the TLL re- 
gion (see FIG.H). The finite size effect is very small, but 
it has an 0{L^^) size dependence. Thisj-eprrection is 
explained by the irrelevant fields {x — 4)cZll23 



L_2i-2l, {(£-2)' + (i-2)'}l. 



(26) 



Although the compressibility should be non-negative in 
the thermodynamic limit (L — > oo), there exists a region 
with K~^ < 0. The reason for this phenomenon is con- 
sidered as follows: we have done the computation under 
the micro canonical ensemble, while the theory of TLL is 
constructed in the grand canonical ensemble. The nega- 
tive inverse compressibility is expected to approach for 
larger systems (see FIG.^. 

FIG. H shows that the Drude weight decreases to zero as 
J/t is increased, but the value remains finite at the phase 
boundary J = Jc determined by k~^ = 0. There is also 
an 0{L^^) size dependence for J < J^. However, for J > 
Jc, the slopes of the Drude weight become steeper as the 
system size increases, and it seems to have a discontinuity 
in the thermodynamic limit. This is consistent with our 
discussion that the Drude weight remains finite at the 
phase-separation boundary. 

In order to check the consistency of the relation for the 
TLL, we compare the charge velocities obtained by the 
two independent methods (see FIGJJ); one is obtained 
from the low-energy spectrum as 



Ei{L;N;S = 0)-EomN) 

27r/L 



(27) 



where Ei{L;N;S = 0) is the first singlet excited state 
for the wave number 27r/L, and the other is given by the 
relation derived from (|l6| ) and ([l7|), 



2D 



(28) 



We have seen that the Drude weight and the compress- 
ibility have the little size dependence for J < Jc. It is 
natural that the size dependence of the charge velocity 
obtained by ( p7| ) is larger than that of (p8|), since the for- 
mer is given by the difference of the energies for discrete 
wave numbers, but the latter is given by the differentia- 
tion of the energy by the continuous quantity <I>, so that 
we used finite size data (at L = 16) for (p8|). We find 
that the charge velocities calculated in the two ways are 
consistent in the region J/t = 2.0-3.0, extrapolating ( |27| ) 
as Vc{L) = vcioo) + Ajl? + B/L^. 

Finally, to see the degeneracy in the phase-separated 
region, we do the Legendre transformation / = e — /xn, 
where e = -Eq/-^- The chemical potential /x is defined by 



KL;N) 



Eo{L;N + 2)-Eo{L;N-2) 



(29) 



The / versus J/t at quarter-filling is shown in FIG.^ As 
the system size L is increased, there appear two regimes 
with different slopes, and the value / approaches in 
the phase-separated region (see also FIG.|). This is the 
evidence that the transition is the first order. In fact, 
the phase-separation boundary Jc can also be defined at 
the point of_£ = which is equivalent to the Maxwell 
constructiorO. In finite systems, the Maxwell construc- 
tion tends to estimate the phase bomidary at smaller J 
than the one determined by k~^ = oIJ'lHI. 

FIG.^ shows that / is almost constant against the 
change of the electron density in the phase-separated re- 
gion. This result suggests that the phase-separated state 
approximately has degeneracy 0{L), which corresponds 
to the large degeneracy 0{L) in the ferromagnetic region 
(A > 1) of the spin S XX Z model. 



IV. SUMMARY AND DISCUSSIONS 

In this paper, we have studied the instability of the 
c — 1 CFT for the U(l) symmetry case. For the spin- 
less fermion (or the S — 1/2 XX Z spin chain), con- 
sidering that the excitation energies are limited by the 
band-width, we showed that the sound velocity is in- 
verse to the Gaussian coupling v ex 1/K. In addition, 
since the transition from the metallic (XY) phase to the 
phase-separated (ferromagnetic) state is the first order, 
the spin wave velocity and the Gaussian coupling behave 
as vK ~ Const., v/K ex 1 — A. 

About the ID t-J model, since the low-energy behav- 
ior of the charge part is described by the c = 1 U(l) 
CFT, and the critical exponent Kp diverges near the 
phase separation, we expected that the compressibility 
is proportional to {Jc — J)~^, and the Drude weight is 
constant. The obtained numerical results supports this 
expectation. 

These asymptotic behaviors of the TLL have been ob- 
served in many cases; the linear behavior of theJuverse 
compressibility was also observed by Troyer et a/.til in the 
ID extended t-J model with the (next-)nearest-neighbor 
repulsion. The discontinuity of the Drude weight (kinetic 
energy) was found by Sandwik and Sudb0 in the case of 
the ID two band Hubbar|d-| model, with increasing the 
nearest-neighbor repulsionE3. Our argument gives a uni- 
fied interpretation of these behaviors. In particular, the 
latter phase-separated state is complicated, |-and similar 
to the one which was found by Sano and OnoQ. Although 
these phase separations are different from the ID t-J 
model, the universality class is expected to be the same, 
so that our argument can be applied to these cases, and 
the discontinuity of the Drude weight and other asymp- 
totic behavior can be explained in the same way. 

Historically, the similar asymptotic behavior in K ^ 
oo was discussed by the g-ologyO. However, in the exactly 
solvable electron systems such as the simple Hubbard 
model or the supersymmetric t-J model {J/t — 2), there 



is no phase separation, therefore such a possibility has 
not been considered seriously. 

Although our speculation is successful in explanation 
of the behavior of the TLL with K ^' oo, the higher 
corrections to vK, v / K near the critical point arc quite 
ambiguous. There are complicated effects come from the 
surface energy. It will be useful to investigate this behav- 
ior by analyzing the t-V model where the exact solution 
is available. 
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{(j)\SI S^ I (j)} = Const . for any i^j, ( A6) 

therefore. 



^z\2| 



*#i 



z\2\ 



L{L-lMS^SLM+LmS^) 



(A7) 



Considering {(f)\{S^)^\(f))/{4>\4>) oc S^, we obtain 
(0|S';S'|+i|(/))/(0|0) oc -S'^/L. Then, the first order per- 
turbation is 

E, = {cj,\ni\<j>)/{m 

= -(A-l)^(0|5|5|+i|<^)/(0|0) 
cx-(A-l)S'2. (A8) 

Therefore, the energy gap between the singlet ground 
state and the fully ferromagnetic state is 

AEo.±sL = S^L{1 - A)[l + Oil - A, 1/L)]. (A9) 



APPENDIX A 



Here we show that the relation (H) applies not only to 
the integrable case S = 1/2 but also to the non-integrable 
case, 5* arbitrary. 

Close to A = 1, we can estimate the ground state 
energy by the perturbation. We treat the SU(2) ferro- 
magnetic term 

no^-Y.S,■S,+^, (Al) 

as a free part, and the anisotropic part 

Hi = -(A-l)^5|5|+i, (A2) 

as a perturbation. 

The energy of the fully ferromagnetic state S^ = ±SL 
is exactly 



E 



ferro 



= -AS'L. 



(A3) 



For A < 1, the ground state has the quantum number 
S^ = 0,q — 0. At A = 1, the ground state wave function 
in the S^ — 0,q — space is derived from the fully 
ferromagnetic state 



{S^)''-\S^^SL), 



(A4) 



where S^ = ^ Sf . The zero-th order energy is given by 

Ho\cl>)=-S^L\cl)). (A5) 

To calculate the first order perturbation, we have to eval- 
uate {(l>\SjSLi\(l>). Since \(j)) is invariant under the per- 
mutation of the lattice sites {j}, we obtain 
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FIG. 1. Compressibility k as a function of J/t for different 
system sizes at n = 1/2. The size dependence at J/t = 2.5 is 
shown in the inset. 
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II L 
FIG. 2. Size dependence of the inverse compressibility in 
the phase-separated regime at n = 1/2, J/t = 4. 




FIG. 3. Drude weight D as a function of J/t for different 
system sizes at n = 1/2. The arrow indicates the critical point 
determined by kT^ = 0. The size dependence at J/t = 2.5 is 
shown in the inset. 
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FIG. 4. Charge velocity Vc derived by (ETl) for different 
system sizes, the extrapolated one (lines), and the one derived 
from (psl) by 16 site data (x), as a function of J/t at n = 1/2. 




FIG. 5. Legendre transformed ground state energy den- 
sity / = e — /xn as a function of J/t for different system 
sizes at n = 1/2. The point / = corresponds to the 
phase-separation boundary determined by the Maxwell con- 
struction. 




FIG. 7. Legendre transformed ground state energy density 
/ = e — /in as a function of the electron density n in L = 16 
system for various values of J/t. 
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FIG. 6. Size dependence of the f = e 
phase-separated regime at n = 1/2, J/t — 4. 
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